On the Quantum Computational Complexity of the Ising Spin Glass Partition 

Function and of Knot Invariants 
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It is sliown that tlie canonical problem of classical statistical thermodynamics, the computation 
of the partition function, is in the case of ±J Ising spin glasses a particular instance of certain 
simple sums known as quadratically signed weight enumerators (QWGTs). On the other hand it 
is known that quantum computing is polynomially equivalent to classical probabilistic computing 
with an oracle for estimating QWGTs. This suggests a connection between the partition function 
estimation problem for spin glasses and quantum computation. This connection extends to knots 
and graph theory via the equivalence of the Kauffman polynomial and the partition function for the 
Potts model. 



I. INTRODUCTION 



Feynman famously conjectured that unlike classical 
computers, quantum computers should be efficient at 
simulating quantum mechanics 1, 2j. This was verified 
by Lloyd 3], and further developed by several other au- 
thors, who demonstrated exponential speedup over the 
best known classical algorithms for a variety of quantum 
mechanical problems such as solution of the constant po- 
tential Schrodinger and Dirac equations using lattice gas 
automata ^,^§^,2], solution of the Schrodinger equation 
in the circuit model simulation of fermionic sys- 

tems ^201 J computation of the thermal rate constant of 
chemical reactions , computation of correlation func- 
tions [l3|, simulation of topological field theories [l5l |. 
and simulation of pairing Hamiltonians such as the BCS 
model fl^ . A naturally related question is whether quan- 
tum computers can efficiently solve problems in classical 
physics. This question was first raised, and partially an- 
swered, in the context of Ising spin glasses 0. It has 
recently received renewed attention in the context of hy- 
drodynamics 0, (with polynomial speedups) , chaos 
|20l I21L 12^ (with exponential speedups, though some of 
these have been challenged [23), and knot theory 
(so far without speedup), which has a deep connection 
to classical statistical mechanics psj . 

Here we revisit classical Ising spin glasses and also ad- 
dress knot theory. The canonical problem of classical sta- 
tistical thermodynamics is the calculation of the partition 
function Z . For a system in thermodynamic equilibrium, 
if the partition function is known, one can obtain ex- 
act results for all thermodynamic quantities such as the 
magnetization, susceptibility and specific heat. Models 
for which analytical calculations of this type can be per- 
formed include a variety of one dimensional (ID) models 
and the 2D Ising model ^^^J. However, for most sys- 
tems of interest, including the 3D Ising model and most 
Ising spin glass models, no analytical calculation of the 
partition function is available [23 ■ We consider classical 
spin systems, in particular the Ising model |^ in which 
each spin has two states cr^ = ±1, and spins interact 



pairwise with an interaction energy of the form JijOiUj. 
From a computational complexity perspective this pro- 
vides a rich class of problems. In particular, the problem 
of finding the ground state of the short range 3D Ising 
spin glass (quenched random J^), as well as the fully 
antiferromagnetic (all Jij = — |J|) 2D Ising model in the 
presence of a constant magnetic field was shown by Bara- 
hona to be NP-hard, by a mapping to problems in graph 
theory. "'^ On the other hand, it is known that there exists 
a fully polynomial randomized approximation scheme for 
the ferromagnetic Ising model (all Jij = | J|), on arbitrary 
graphs ■ 



The problem of sampling from the Gibbs distribution 
of the Jij — ±J (with random signs) spin glass on a 
quantum computer was addressed in |17l |. A linear-time 
algorithm was found for the construction of the Gibbs dis- 
tribution of configurations in the Ising model, including 
partially frustrated models. A magnetic field can be in- 
corporated as well without increase in the run-time. The 
algorithm was designed so that each run provides one 
configuration with a quantum probability equal to the 
corresponding thermodynamic weight. I.e., the probabil- 
ities of measuring states are ordered by the energies of 
the corresponding spin configurations, with the ground 
state having the highest probability. Therefore the par- 
tition function is approximated efficiently and statisti- 
cal averages needed for calculations of thermodynamic 
quantities obtained from the partition function, are ap- 
proximated in the fastest converging order in the num- 
ber of measurements. Unlike Monte Carlo simulations 
on a classical computer, consecutive measurements on a 
quantum computer can be totally uncorrelated. Thus 
the algorithm neither suffers from critical slowing down 
(a polynomial slowdown in Monte Carlo moves associ- 
ated with large correlated spin clusters forming near the 



^ NP-hard problems are those whose proposed solution cannot 
even be verified using a nondeterministic Turing machine in poly- 
nomial time; e.g., is the proposed ground state truly the ground 
state? 
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critical temperature) |3j| , nor gets stuck in local minima. 
This uniform performance is an advantage over the best 
known classical algorithms, which are tailored to specific 
lattices or graphs ^V\. However, the main problem of the 
algorithm is the limited control it offers in the construc- 
tion of a specific realization of bonds on the Ising lattice. 
Indeed, since the run-time of the algorithm is linear, it 
is reasonable to suspect that it cannot simulate a hard 
instance of an Ising spin glass. 

A completely different approach to sampling from the 
Gibbs distribution for the ferromagnetic Ising model was 
recently developed by Nayak, Schulman and Vazirani 
(NSV) [s^, which however does not appear to provide 
a speedup over the best known classical algorithm [33l| . 
The VNS algorithm for the ferromagnetic Ising model 
uses an interesting representation of the partition func- 
tion as a discrete Fourier transform, in conjunction with 
a Markov chain sampling procedure. Two of the main 
results of the present paper are (i) the generalization of 
this Fourier transform representation to the ± J spin glass 
case, (ii) (the central result) a connection of this repre- 
sentation to certain simple sums known as "quadratically 
signed weight enumerators" (QWGTs). Let us now mo- 
tivate the importance of result (ii) . 

In virtually all previous work on simulation of physi- 
cal systems on quantum computers P, II , IS B, la I p , 
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the approach pursued was one of attempting to find a 
concrete algorithm for a specific simulation problem. A 
fruitful alternative is to consider instead the question of 
the complexity class that the simulation problem belongs 
to. We do this here by following a lead due to the Knill 
and Laflamme (KL): KL showed that quantum comput- 
ing is polynomially equivalent to classical probabilistic 
computing with an oracle for estimating QWGTs 34]. 
Combined with our results (i),(ii), this suggests that the 
quantum computational complexity of sampling from the 
Gibbs distribution for the Ising spin glass problem can 
be understood in terms of QWGTs. However, we have 
unfortunately not yet been able to establish the connec- 
tion at this level. Nevertheless, the possibility of a (spin 
glass)-(QWGTs)-(quantum computation) connection is 
sufficiently tantalizing to point it out in detail. We hope 
that by expressing the partition function as a QWGT 
we have taken the first step in a direction that will al- 
low future research to explore the important question of 
the quantum computational complexity of the Ising spin 
glass problem. 

In fact, the connections do not end here. There is a 
rich inter-relation between classical statistical mechanics 
and topology, in particular the theory of classification 
of knots. The first such connection was established by 
Jones 35], who discovered knot invariants (the Jones' 
polynomial) du ring his investigation of topological prop- 
erties of braids |33 . It is known that classical evaluation 
of the Jones' polynomial is #P-hard |33]- The connec- 
tion between knots and models of classical statistical me- 
chanics was greatly embellished by Kauffman 25]. Here 



we will exploit this connection to show [our main result 
(iii)] that the evaluation of another knot invariant, the 
Kauffman polynomial, can also be cast in some cases as a 
QWGT evaluation problem. Thus a quantum algorithm 
for QWGT evaluation should shed light on the quantum 
computational complexity of knot invariants, a subject 
which has been explored by Freedman et al. [3^ and 
by Kauffman and Lomonaco 39]. Knot invariants are, 
in turn, also tightly related to graph theory; e.g., the 
graph coloring problem can be considered an instance of 
evaluation of the Kauffman polynomial, via the Tutte 
polynomial 25] . 

Mathematically, the reason that these seemingly unre- 
lated subjects are all inter-related is due to the fact that 
key properties can be expressed, in all cases, in terms of 
certain polynomials. While these polynomials originate 
from widely distinct problems, from the point of view of 
computational complexity their evaluation is one and the 
same problem, much in the same spirit as the fact that 
solving one problem in the class of NP-complete prob- 
lems solves them all licj. The present work contributes 
to this unification. 

The structure of this paper is as follows. In Sec. ^ 
we derive our first main result: we rewrite the Ising spin 
glass partition function for an arbitrary graph as a dis- 
crete Fourier transform. This motivates the considera- 
tion of the partition function evaluation problem in terms 
of its computational complexity, which we formalize in 
Sec. Uni We then review QWGTs in Sec. |TV| In Sec- 
tion we derive our second main result: the connection 
between the evaluation of the Ising spin glass partition 
function and QWGTs. In Sec. I VII we continue the pro- 
gram of connecting disparate objects to the problem of 
QWGT evaluation, and obtain our third main result: we 
show that the Kauffman polynomial too can be expressed 
a QWGT. We do this after first reviewing the connec- 
tion between knots and classical statistical mechanics. 
Section IVIII concludes. Some further observations con- 
cerning the representation of the partition function are 
collected in Appendix 1X1 

II. FOURIER TRANSFORM 
REPRESENTATION OF THE PARTITION 
FUNCTION 

Let G = {E, V) be a finite, arbitrary undirected graph 
with \E\ edges and \V\ vertices. Identify each vertex 
i € V with a classical spin {ai = ±1) and each edge 
(j,j) G E with a bond {Jij — ±J). Denote a given 
spin configuration by cr = [ai, (72, f^lv]) ^^^d a bond 
configuration by (J12, Jy, ...). We assume that the 
bond configuration is chosen at random and then remains 
fixed ( "quenched randomness" ) . The Hamiltonian of the 
system is 

F(cr) = - ^ JijCTiCTj- (1) 
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(We remark on the case with a magnetic field in Ap- 
pendix The probabihty of the spin configuration a 
in thermal equilibrium at temperature T is given by the 
Gibbs distribution: 



P{<y) = -^W{a), (2) 

where the Boltzmann weight is 

W{<j)^eM^PH{cj)], (3) 

(3 = l/kT is the inverse temperature, and Z is the parti- 
tion function: 

Z{{J.,})=J2eM-PH{a)]. (4) 
Now note the identity 

= cosh(a;)[l -t- tanh(.T)] (5) 
and use it to rewrite the Boltzmann weight Q as 

W{a)^ n cosh(/3J,jCT,crj)[l +tanh(/3jycr,;crj)]. (6) 

Let 

where qij — ±1 is a quenched random variable. Since 
cosh(a;) = cosh(— x) and tanh(a;) = — tanh(— a;) we find 

W{a)=Q W (l + g,,a,CT,A), (8) 

where 

e= [cosh(/3J)]l^l, (9) 

and 

A = tanh(/3J). (10) 
Next expand out the product to obtain 

P{a) = e[l + A (JiqijCFj 

(■i,i)eB 

+A^ ^ (ai(7jjCrj)((Tfc(7fcicri) 4- ....](11) 

Note that A*"' is the coefficient in front of a sum containing 
k bonds qij, which are not necessarily all connected. For 
example, the term qijCJiCJjqkiUkUi where all indexes differ 




(a) (b) (c) 



FIG. 1: Illustration of subgraphs. 

represents two unconnected bonds (Fig. whereas if 
j = k the same term represents two connected bonds 
sharing one spin (Fig. ^). Let b C E denote such a 
subgraph, with k = \b\ edges. There are ('|^|') ways of 
choosing subgraphs with |6| edges. Thus the total number 
of subgraphs is: 




This suggests to label the subgraphs in a binary fash- 
ion (to be explained below): Let b = (6i, 62, 6|£;|) 
(where bj G {0,1}) be the binary number of subgraph 
b. The numbering is such that bj = 1 (0) indicates the 
presence (absence) of edge number j of G. Further, we 
will use the convention that first all \E\ single-edge sub- 
graphs are counted (i.e., vectors b with a single 1 en- 
try, all the rest 0), then all ('^') double-edge subgraphs 
(vectors b with two 1 entries, all the rest 0), etc. Thus, 
e.g., b = (0,0, ...,0, 1) could corresponds to the subgraph 
containing only the first edge (or bond): J12, whereas 
b = (0, 0, ...0, 1, 0, 1) could correspond to the subgraph 
containing only J12 and J34. Note that |6| is the Ham- 
ming weight of b. Since the total number of subgraphs is 
2l^l, the above numbering scheme is a one-to-one cover- 
ing of the subgraph space, and moreover, the subgraphs 
are labeled in increasing order of 

Next note that in Eq. Ijlll) , a spin at may appear more 
than once in sums of order A; > 2, in fact as many times 
as the number of bonds emanating from ct^, which we 
denote by deg^(i). Then cr^^StC*) ^^iq contribution of 
spin i in subgraph b to the product in the sum of order 
k ^ |6| in Eq. (|ll|l . Collecting all the observations above, 
it follows that Eq. ffTT^ can be rewritten as: 

p{a)=eY^x^'^ n ^^'^'^ n (13) 

b ieV{b) {i.j)eb 

Here V{b) denotes the set of vertices in subgraph b. 
Next introduce the parity vector a'' = {a\, a|'y|) 

of a subgraph, with components = 1 if (i G V{b) and 
degt(z) is odd), = otherwise. Then clearly (7^°Sb(«) _ 

tj"' . At this point it is more convenient to transform to 
a binary representation for the spins as well. Let: 
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be the components of the binary vector of spin values s 
(si,S2,...,S|y|). We have = (-1)^' so that crf*^fW 
Therefore: 



n - 

ieVib) 



dogb(i) 



(-1)" 



(15) 



where "■" stands for the (mod 2) bit-wise scalar prod- 
uct. The same change can be affected for the bonds by 
introducing: 



Wi 



1^3 



(16) 



so that the binary vector w = (1012,^13,...) of length 
\E\ specifies whether edge supports a ferromagnetic 
(wij — 0) or antiferromagnetic {wij — 1) bond. Again, 
= (-1)™'J , so that: 



n *J 



(-1) 



b-w 



(17) 



Using Eqs. 1^511,(271) in Eq. it is now possible to 

rewrite the Boltzmann weight of a particular spin config- 
uration s as: 



W{s) ^ e^Ai''i(-i)"'- 



s-\-b-'w 



(18) 



Now, the partition function is just the sum over all spin 
configurations. Using Eq. (|18() we thus find: 

Proposition 1 The spin-glass partition function is a 
double discrete Fourier (or Walsh-Hadamard) transform, 
over the spin and subgraph variables; 



■s-\-b-w 



(19) 



s b 



Note that the sum over s extends over 2l^l terms, while 
the sum over b extends over 2l^l terms. By changing 
the order of summation, the sum over s can actually be 
carried out, since: 



i-ir 



(20) 



where 6 is the Kronecker symbol. A systematic procedure 
for finding the parity vectors from the subgraphs vectors 
uses the incidence matrix A. For any graph G = {E, V) 
this matrix is defined as follows {v G V) 41.] : 



A _ , 1 (w = i and (i,j) G E) , . 

^v,{^,J) - \ else ' ^^^> 



so A is a |y| X \E\ matrix of O's and I's. It is well known 
that given G and a specific subgraph b |4ll |. 



a" ^A-b. (22) 
Combining Eqs. @ , ((101) , CH , (EOJ , (EH) , we finally have: 



Z{w) 



2l^l 



(l_^2)|£|/2 



b-w 



(23) 



6: A-b=0 



In words, the sum over the subgraphs includes only those 
with zero overall parity, i.e., those having an even num- 
ber of bonds emanating from all spins. This immediately 
implies that "dangling-bond" subgraphs are not included 
in the sum. We note that Z(w) can also be rewritten as a 
power series in A, which is useful for a high-temperature 
expansion; this is discussed in App. lA2l The representa- 
tion (|23() allows us to establish a direct connection with 
QWGTs, which is the subject of the quantum computa- 
tional complexity of Z, to which we turn next. 



III. FORMULATING THE COMPUTATIONAL 
COMPLEXITY OF THE ISING SPIN GLASS 
PARTITION FUNCTION 

The most natural computational complexity class 
for quantum computation is BQP: the class of deci- 
sion problems solvable in polynomial time using quan- 
tum resources (a quantum Turing machine, or, equiv- 
alently, a uniform family of polynomial-size quantum 
circuits) with bounder probability of error [4^ |4^ . 
The class BQPP is the natural generalization of BQP 
to promise problems [s^l- Relative to the polynomial 
hierarchy of classical computation, it is known that 
BPPCBQPCPPCPSPACE, but none of these inclusions 
is known to beproper |4^ . 

In order to address the quantum computational com- 
plexity of the spin glass partition function we define: 

Definition 1 An instance of the Ising spin glass problem 
is the data A = (G, w, A). 

Now, let Z^{s) be the partition function for given data 
A and given spin configuration s G {0, 1}'^', and let 



Zf(s):{0,ir-{0,1} 



(24) 



be the /c**^ digit of Z^{s). Let r^(s,fc) be a quantum 
circuit that takes as input the spin configuration s and 
the digit location fc, for fixed data A. Let Wp — Yii Ui be 
an implementation of P'^ in terms of some unitary oper- 
ators Ui. Let the circuit be designed so that the answer 
Z^[s) is encoded into the state of the first qubit, and let 
III = |0)i(0| be the projection onto state |0) of this qubit. 
Then the probability of measuring the state |0)i after the 
circuit was executed, starting from the "blank" initial 
state |0) = I0---0), is Pr[r^(s,fc)] = (OlVFjJniWrlO). 
We can now define: 
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Definition 2 Z^{s) S BQP if there exists a classical 
polynomial-time algorithm for specifying such that 

Pr[r^(s,fc)] > I if Zt{s)=Q and 

Pr[r^(.,fc)] < i if Zt{s) = l. 

We can then formulate the foUowing open problem: 

Problem 1 For which instances A of the Ising spin glass 
problem is evaluating the partition function in BQP? 

A particularly promising way to attack this problem 
appears to be the connection to QWGTs, which we ad- 
dress next. 



IV. QUADRATICALLY SIGNED WEIGHT 
ENUMERATORS 

Quadratically signed weight enumerators (QWGTs) 
were introduced by Knill and Laflamme in Ref. |3J| 
(where "QWGT" is pronounced "queue-widget"). A gen- 
eral QWGT is of the form 

S{A,B,x,y)^ ^ (-l)''^^^!''!?;"-!''! (25) 

b:Ab=0 

where A and B are 0, 1-matrices with B of dimension 
n X n and A of dimension m x n. The variable b in the 
summand ranges over 0, 1-column vectors of dimension n, 
and all calculations involving A, B and b are done modulo 
2. It should be noted that \S{A,B,x,y)\ < {\x\ + |y|)". 
In Ref. [3^ it was shown that quantum computation is 
polynomially equivalent to classical probabilistic compu- 
tation with an oracle for estimating the value of certain 
QWGTs with x and y rational numbers. In other words, 
if these sums could be evaluated, one could use them to 
generate the quantum statistics needed to simulate the 
desired quantum system. 

More specifically, let / be the identity matrix, diag(74) 
the diagonal matrix whose diagonal is the same as that of 
A, and Itr(A) a matrix formed from the lower triangular 
elements of A (the matrix obtained from A by setting to 
zero all the elements on or above the diagonal). Then 
for: 

Problem 2 KL promise problem; Determine the sign of 
S(A,\tr{A),k,l) with the restrictions of having A square, 
diag(^) = I, k and I being positive integers, and the 
promise \S{A,\tT:{A),k,l)\ > {k^ + P^/'^ /2. 

KL demonstrated the following: 

Theorem 1 (Corollary 12 in ^]): The KL promise 
problem is BQPP-complete. 

KL's strategy in showing the connection between 
QWGT evaluation and quantum computation was to 
show that in general expectation values of quantum cir- 
cuits can be written as QWGTs. 



V. THE PARTITION FUNCTION - QWGT 
CONNECTION 

We are now ready to prove our central result. 

Theorem 2 The spin-glass partition function is a spe- 
cial case of QWGTs. Specifically: 

^^1^;^ Z{w) = S{A,dg{w),X,l). (26) 

Here dg{w) is the matrix formed by putting w on the di- 
agonal and zeroes everywhere else, and A is the incidence 
matrix of G. 

Proof. In Eq. (|25|l identify b as the subgraphs of G = 

(E, V), n = \V\, m = \E\, and note that when B = dg(u') 

bBb = biWibi — b ■ w 

i 

since bi = or 1. Then Eq. H26|l follows by inspection of 
Eqs. and (gni). ■ 

Corollary 1 Evaluating the spin glass partition function 
is in #P? 

Proof. The problem of evaluating QWGTs at integers 
is in the class #P [s^].^ In our case x = 1, y = X and the 
coupling constant J can always be chosen so that A = 
tanh(/3J) is integer. ■ 

It is tempting to check the relation of Theorem 13 to 
the KL promise problem (Problem [21 . It follows from 
Eq. ijini), from Z > 0, and from < A = tanh(/3J) < 1, 
that sign[S'(yl, dg(it;). A, 1)] = Hence, unfortunately, 
the KL problem in its present form is of no use to us. 

Further consideration reveals that, while the constraint 
that k and / are positive integers is easily satisfied, 
and the promise takes a nice symmetric form: Z{w) > 

2'^'~"'^-|y3^2^W72; the remaining constraints - A square, 
diag(y4) — I, B — Itr(^) ~ anyhow result in severely re- 
stricted instances of spin glass graphs. We thus leave as 
open the following problem, inspired by the KL problem: 

Problem 3 Formulate a promise problem in terms 
of Z{w) [or, equivalently, S{A,dg(w),X,l)] which is 
BQPP-complete. 

We turn next to showing the connection between our 
discussion so far and problems in knot theory. 



^ #P is the class of function problems of the form "compute f{X)" , 
where / is the number of accepting paths of a nondeterministic 
polynomial-time Turing machine. The canonical #P-complete 
problem is #SAT: given a Bo olean f ormula, compute how many 
satisfying assignments it has l4(l 1431 . 

^ It contains the problem of evaluating the weight enumerators of 
a b inary linear code at rational numbers, which is #P-complete 
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FIG. 2: Definition of crossing variable for knots. 



VI. THE PARTITION FUNCTION - KNOTS 
CONNECTION 

The canonical problem of knot theory is to determine 
whether two given knots are topologically equivalent. 
More precisely, in knot theory one seeks to construct a 
topological invariant which is independent of the knot 
shape, i.e., is invariant with respect to the Reidemeister 
moves |25| . This quest led to the discovery of a num- 
ber of "knot polynomials" (e.g., the Jones and Kauffman 
polynomials) 25]. These also play a major role in graph 
theory as instances or relatives of the dichromatic and 
Tutte polynomials e.g., in the graph coloring prob- 
lem. Roughly, two knots are topologically equivalent iff 
they have the same knot polynomial. It is well known 
[25L Issf that there is a connection between knot polyno- 
mials and the partition function of the Potts spin glass 
model, a generalization of the Ising spin glass model to 
q > 2 states per spin: 

Hpottsiq, s) = - ^ JijSsi,sj, (27) 

where Si G {0,...,q — 1} and (5s.^s^ — 1 (0) if s,; = sj 
{si ^ Sj). We first briefly review this connection. 

Consider a knot embedded in 3D-space (imagine, e.g., 
a piece of rope). In the standard treatment j2^, the knot 
is projected onto the plane and one obtains a "2D-knot 
diagram" . The essential topological information about 
the knot is contained in the pattern of "crossings", the 
2D image of where one rope segment went over or under 
another rope segment. A crossing takes values 6^ = ±1 
according to Fig.|21 47]. 

A connection to spin glasses can be made by assign- 
ing quenched random values to the crossing variables bk, 
so that the links cross above and below at random. It 
was shown by Nechaev that in this case the Kauff- 
man polynomial is identical, up to an irrelevant multi- 
plicative factor, to the Potts model partition function, 
Zpottsiq A Jij}) = Z]sexp[-/3ffpotts(g, s)]- To explain 
this connection we need to introduce some terminology. 
The 2D-knot diagram lives on a lattice M composed of 
lines oriented at ±45°, intersecting at the crossings bk, 
that carry the disorder. One can define a dual lattice C, 
rotated by 45°, so that its horizontal and vertical edges 
(denoted bij) are in one-to-one correspondence with the 
vertices bk of 7W (Fig. 6 in Q). Let 

^ ( — 6fe if the (ij)-edge is vertical , ^ 

1 bk if the (ij)-edge is horizontal ' ^ ^ 



The Potts spin states Si are connected to knot properties 
in an abstract manner; they are related to the Kauff- 
man polynomial variable A, which in turn is a weight for 
the manner in which 2D-knot diagram is disassembled 
into a set of microstates (and also related to the Jones' 
polynomial variable t: the Jones and Kauffman polyno- 
mials coincide when t = A^^'^). Precise definitions can 
be found in for our purposes what matters is that 
the equivalence of the Kauffman polynomial to the Potts 
spin glass partition function is established once one as- 
signs the Potts variables q and f3Jij the values 

q = {A^+A-y, 
(3 J,, = ln(-A-4f..). (29) 

With these identifications Nechaev has shown that 
the Kauffman polynomial (knot invariant) {K{A)) ~ 
c{A, {bij})Zpottsi,q, J), where the constant c does not de- 
pend on the spin states Solving for A we find 
A = ±[iq^/^ ±iq- 4)i/2)/2]i/2. Thus /3Jy can be real 
only for g > 4. In the Ising spin glass case {q — 2) 
we obtain a complex- valued f3Jij, which in turn implies 
complex- valued A = tanh(/3J), and hence the estimation 
of the QWGT polynomial with complex- valued x. 

Finally, we note that a physically somewhat unsatis- 
factory aspect of the knots-Potts connection is that now 
the (complex-valued) temperature cannot be tuned inde- 
pendently from the bonds Jij. However, this docs not 
matter from the computational complexity perspective: 
we have established our third main result: 

Proposition 2 Computing the Kauffman polynomial at 
q = 2 is equivalent to the problem of computing the 
QWGT polynomial with complex-valued x.'^ 

Hence an efficient quantum algorithm for estimating 
QWGTs will be decisive for knot and graph theory as 
well. 



VII. CONCLUSIONS 

The connection between QWGTs and quantum com- 
putational complexity established by KL on the one 
hand, and the connection between QWGTs and the spin 
glass and knots problems established here on the other 
hand, suggests that the quantum computational com- 
plexity of spin glass and knots problems may be decided 
via the connection to QWGTs. Similar remarks apply 
to a number of combinatorial problems in graph theory, 
via their well-established connections to knot theory. In 
particular, it would be desirable to find out the quantum 



* This knots-QWGT connection does not appear to hold in the 
case q > 2, since in this case we cannot separate Ss^,sj into a 
product of single-spin variables, a step that is essential in deriving 
the representation of Z as a QWGT [see text around Eq. 1131 ]. 
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computational complexity of questions framed in terms 
of properties of S{A, dg(w), A, 1), with A real (Ising spin 
glass) or complex (Kauffman polynomial with q = 2). 
We leave these as open problems for future research. 
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2. Power Series Representation 



Another useful representation of Eq. H23|) can be ob- 
tained by grouping together all subgraphs with the same 
number of edges. To this end, let fej*^-* denote the j*'^ sub- 
graph with k edges. According to the numbering scheme 
introduced in Sec. m the corresponding binary number 
of such a subgraph, bj'^^ is the j*"^ permutation of a vec- 
tor of exactly k I's and \E\ — k zeroes. Since |6| is the 
Hamming weight of b these subgraphs all have |6| — k. 
There are ('^') such subgraphs, all with equal weight A'^. 
Therefore: 



APPENDIX A: ADDITIONAL OBSERVATIONS 

1. The case with a Magnetic Field 

A magnetic field B can be included in the Hamiltonian 
[Eq. Ql], by adding a term — X^iey Biai. We can repeat 
the analysis above by introducing a fictitious "always- 
up" spin, numbered 0. In this manner we can rewrite the 
magnetic field term as 



J2 B^a^ - ^ S 



(Al) 



where = \. The corresponding graph has a "star" 
geometry, with spin in the center, connected to all other 
spins, which in turn arc connected only to spin (Fig.^). 
The analysis above can then be repeated step- by-step, 
with the relevant subgraphs being those of the star graph. 
However, we then cannot recover the QWGT form, due 
to the extra summation over the star-graph subgraphs: 
Denote the latter b' . Since each spin in V is connected 
once to the central spin ctq, the star graph subgraphs all 
have trivial parity vectors, a\ =1. Then Eq. H2U|I is 
replaced by 



(q''+q'' ,0)- 



(A2) 



This causes a violation of the condition h : Ab = {) needed 
in the definition of the QWGT sum. Thus it appears that 
QWGTs do not include the case with a magnetic field. 



I 



(l_X^)\E\li 



h. -be 

d (fc) 

(q J ,0) 



(A3) 



In this form we have a series expansion in powers of A, 
corresponding to the number of edges of the subgraphs. 

A clear simplification results in the fully ferromagnetic 
Ising model (^+), where w = (0, 0, 0), and in the fully 
antiferromagnetic case {Z-), where w = (1,1,..,1). In 
the latter case we have simply b ■ w = so that com- 
bining the two cases we obtain from Eq. ^Z'6\ : 



Z± = 



(1-A2)|s|/' 



E(±^)'"v 



.0)- 



(A4) 



Eq. (jASj), on the other hand yields: 



Z+ - 



2l^l 



(l-A2)|s|/2 



i+E(±^)'E^, 



fc=l 



(fc) 

( a 3 ,0) 



(A5) 

As already remarked, there exist an efficient classical al- 
gorithm for calculating Z in the case of the fully ferro- 
magnetic Ising model 
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